Conduction eigenchannels of atomic-sized contacts: 
Ab initio KKR Green's function formalism. 
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We develop a formalism for the evaluation of conduction eigenchannels of atomic-sized contacts 
from first-principles. The multiple scattering Korringa-Kohn-Rostoker (KKR) Green's function 
method is combined with the Kubo linear response theory. Solutions of the eigenvalue problem for 
the transmission matrix are proven to be identical to eigenchannels introduced by Landauer and 
Biittiker. Applications of the method are presented by studying ballistic electron transport through 
Cu, Pd, Ni and Co single-atom contacts. We show in detail how the eigenchannels are classified 
in terms of irreducible representations of the symmetry group of the system as well as by orbital 
contributions when the channels wave functions are projected on the contact atom. 

PACS numbers: 73.63.Rt, 73. 23. Ad, 75.47.Jn, 73.40.Cg 
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I. INTRODUCTION 

The invention of the scanning tunneling microscope 1 - in 
1981 and a consequent development in the beginning of 
the nineties of the remarkably simple experimental tech- 
nique known as mechanically controllable break junction 
(MCB J)2i2 led to the possibility of fabrication of metallic 
point contacts approaching the atomic scale. The recent 
review article (Ref. |J) summarizes the numerous achieve- 
ments in this field. In the experiments the conductance 
measured as a function of the elongation of the nanocon- 
tacts decreases in a stepwise fashion 3 -^! 7 - with steps of 
order of the conductance quantum Go = 2e 2 /h. Such 
behavior of the conductance is attributed to atomic rear- 
rangements that entails a discrete variation of the contact 
diameter 

The electron transport in metallic nanocontacts is 
purely ballistic and phase-coherent because their size is 
much smaller than all scattering lengths of the system. 
According to Landauer ji 2 . conductance is understood as 
transport through nonmixing channels, 

2e 2 A 

n=l 

where T„'s are transmission probabilities. They are de- 
fined as eigenvalues of the transmission matrix tt^ . Here 
the matrix element r nm gives the probability amplitude 
for an incoming electron wave in the transverse mode 
(channel) n on the left from the contact to be transmitted 
to the outgoing wave in the mode m on the right. Conse- 
quently, the eigenvectors of tt^ are usually called eigen- 
channels. It was shown in the pioneering work by Scheer 
et alJ£ that a study of the current- voltage relation for the 
superconducting atomic-sized contacts allowed to obtain 
transmission probabilities T ra 's for particular atomic con- 
figurations realized in MCB J experiments. The T„'s are 
found by fitting theoretical and experimental I —V curve 
which has a peculiar nonlinear behavior for supercon- 



ducting contacts at voltages eV smaller than the energy 
gap 2 A of a superconductor—. The origin of such effect is 
explained in terms of multiple Andreev reflections p« The 
analysis of MCBJ experiments within the tight-binding 
(TB) model suggested by Cuevas et aZ i 15 ' 16 gave a strong 
evidence to the relation between the number of conduct- 
ing modes and the number of valence orbitals of a contact 
atom. 

To describe the electronic and transport properties of 
nanocontacts, quite a big number of different methods 
which supplemented each other were developed during 
the last 15 years. Early models employed a free-electron- 
like approximation . 9 ' 17 : 18 Further approaches based on 
density functional theory (DFT) used psuedopotentials 
to describe atomic chains suspended between jellium 
electrodesi 19,20 The TB models were applied to the prob- 
lem of the conduction eigenchannel o 15 : 21 and to the study 
of the breaking processes of nanowires.— The up-to-date 
fully self-consistent ab initio method a 23 ' 24 : 25 allowed to 
treat both the leads and the constriction region on the 
same footing and to evaluate the non-equilibrium trans- 
port properties as well 2 - ' 25 : 26 . 

The scattering waves, underlying a concept of eigen- 
channels introduced by Landauer and Biittiker^ do not 
form an appropriate basis for the most of ab initio meth- 
ods. Instead, one considers conduction channels as eigen- 
vectors of some hermitian transmission matrix written in 
terms of local, atom centered basis seti 15 : 24 One of the 
goals of the present paper is to establish a missing link 
between these approaches. Below we introduce a formal- 
ism for the evaluation of conduction eigenchannels, which 
combines an ab initio Korringa-Kohn-Rostoker (KKR) 
Green's function method 27 for the electronic structure 
calculations and the Baranger and Stone formulation of 
the ballistic transport 28 . In recent publications , 29 : 30 we 
have successfully applied this method to the study of the 
electron transport through atomic contacts contaminated 
by impurities. In the present paper, mathematical as- 
pects of the problem are considered, followed by some 
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applications. In particular, we analyze the symmetry of 
channels and relate our approach to the orbital classifi- 
cation of eigenmodes introduced by Cuevas et o/.— 

The paper is organized as follows. A short descrip- 
tion of the KKR method is given in Sec. II. We pro- 
ceed in Sec. Ill with a formal definition of eigenchan- 
nels for the case of realistic crystalline leads attached to 
atomic constriction. Sec. IV supplemented by Appen- 
dices A and B contains mathematical formulation of the 
method. Briefly, using the equivalence of the Kubo and 
Landauer approaches for the conductance , 28 ' 31 we build 
the transmission matrix rr' in the scattering wave repre- 
sentation. The angular momentum expansion of the scat- 
tering Bloch states within each cell is used further to find 
an equivalent, KKR representation of the transmission 
operator for which the eigenvalue problem can be solved. 
Applications of the method are presented in Sec. V. In 
particular, we focus on transition metal contacts (such as 
Ni, Co and Pd), since experimenta l 32 ' 33 ' 34 and theoreti- 
cal studie s 35 ' 36 ' 37 ' 38 ' 39 of their transport properties have 
been attracting much attention during the last years. 
Experiment o 40 ' 41 ' 42 ' 43 ' 44 regarding ballistic magnetoresis- 
tance (BMR) effect in ferromagnetic contacts are com- 
mented. A summary of our results is given in Sec. VI. 



II. ELECTRONIC STRUCTURE 
CALCULATION OF THE ATOMIC CONTACTS 

The systems under consideration consist of two semi- 
infinite crystalline leads, left (L) and right (R), coupled 
through a cluster of atoms which models an atomic con- 
striction. In Fig. 1 a typical configuration used in the 
calculations is shown — the two fee (001) pyramids at- 
tached to the electrodes are joined via the vertex atoms. 
We employed the ab initio screened KKR Green's func- 
tion method to calculate the electronic structure of the 
systems. Since details of the approach can be found 
elsewhere^ 7 - only a brief description is given below. 

In the KKR formalism one divides the whole space into 
non-overlapping, space-filling cells, with the atoms (and 
empty spheres) positioned at the sites R„, so that the 
crystal potential V is expressed in each cell as V n (r) — 
V(R„ + r). The one-electron retarded Green's function is 
expressed in terms of local functions centered at sites R„: 

G+(R„+ r,R„, + r';E) 
= 5 nn ,VEY,Rl(r<;E)H2(r>;E) (1) 

L 

LL' 

where r, r' are restricted to the cells n and n'; r<, r> 
denote one of the two vectors r or r' with the smaller or 
the larger absolute value, and local functions R2(r;E) 
and _ff£ (r; E) are the regular and irregular solutions of 
the Schrodinger equation for the single potential V n (r). 




FIG. 1: Geometry of an atomic constriction: two fee (001) 
pyramids are attached via apex atoms. Conductance is cal- 
culated between the two planes Sl and Sn positioned in the 
leads. 



Here the index L = (I, m) stands for the angular mo- 
mentum quantum numbers and atomic units are used: 
e = — H = 1, m = 1/2. The structural Green's func- 
tion GYL'(E) (structure constants) in Eq. |T|) is related 
to the known structure constants of the appropriately 
chosen reference system by the algebraic Dyson equa- 
tion which includes the difference At = SnnidLL'^l 
between local i-matrices of the physical and a reference 
system. In the screened KKR method 27 we use a lat- 
tice of strongly repulsive, constant muffin-tin potentials 
(typically, ~ 4Ry height) as reference system that leads 
to structure constants which decay exponentially in real 
space. 

Within the screened KKR method both a constriction 
region and the leads are treated on the same footing. 
This is achieved by using the hierarchy of Green's func- 
tions connected by a Dyson equation, so that we per- 
form the self-consistent electronic structure calculations 
of complicated systems in a step-like manner. First, us- 
ing the concept of principal layers together with the dec- 
imation technique^ we calculate the structural Green's 
function of the auxiliary system consisting of semi- infinite 
leads separated by a vacuum barrier. At the second step, 
the self-consistent solution of the impurity problem is 
found by embedding a cluster with perturbed potentials 
caused by the atomic contact into the auxiliary system. 
Due to effective screening of perturbation, the algebraic 
Dyson equation for the structure constants is solved in 
real spaced 

III. DEFINITION OF EIGENCHANNELS 

The concept of eigenchannels is introduced in the Lan- 
dauer approach to ballistic transport, where the problem 
of the conductance evaluation is considered from view- 
point of scattering theory. Following Landauer^ we look 
at the system shown in Fig. 1 as consisting of two semi- 
infinite leads (electrodes) attached to a scattering region 
(atomic-sized constriction). Far away from the scatter- 
ing region the propagating states are the unperturbed 
Bloch waves ^(r, E) of the left (L, z — * — oo) and right 
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(R, z — > +cjo) leads, where k belongs to the isoener- 
getic surface E = const (Fermi surface, E — Ep, in case 
of conductance) and a common notation k — (k, A) is 
used to denote Bloch vector k and band index A. For 
the eigenchannel problem one considers in-coming and 
out-going states in the L and R leads normalized to a 
unit flux. The in-states in L and out-states in R are 
^"fc = ^fc/V^fc with positive velocity Vk oc vt > along 
2-axis. The conjugated states = ^ k * / y/\v_ k \ are 
the out-states in L and in-states in R with negative ve- 
locity V—k oc v z _ k < 0. Here v k is a z-component of 
the group velocity = dEk / dk; a proportionality fac- 
tor between Vk and v k related to a particular choice of 
normalization of the Bloch waves is introduced further in 
Sec.IV.B. 

The potential AV(r) describing the constriction in- 
troduces a perturbation to the perfect conductor. Let 
$fc(r, E) be a perturbed state which is a solution of 
the Lippmann-Schwinger equation for an in-coming state 
in L: 



$ k (r,E) = <t>°k(*,E) 



(2) 



J d 3 r' G+(r,r',E)AV(r')<f>k(r',E) 



where the integral goes over all space, and G^(r, r', E) is 
the retarded Green's function of the perfect conductor. 
Asymptotic behavior of 3>fc(r, E) is 



*k(r,£) 



<Z>° k (r,E) + J2pkk'(E) K*(r,E) 

k' 



k' 



(3) 



where r fcfc ' (E) and p kk > (E) are transmission and reflec- 
tion amplitudes, assuming elastic scattering (Ek = E — 
Ek>). According to the Landauer-Biittiker formula, 12 
conductance is given by g — goTr(rr^), where trace 
goes over in-coming states (fe) in the left electrode and 
E = Ep. An equivalent formulation with respect to in- 
coming states (— k') from the right electrode reads as 
9 = ffoTr^tr). 

Eigenchanncls appear from a unitary transformation 
of in- and out-states. Let w be a unitary transform of in- 
states in L: $° in (r, E)]^^ = Y,k^k{E)^{r,E). 
The corresponding solution 3>„(r) of Eq. @ for an arbi- 
trary r is 



The matrix r, however, is not diagonal in basis v. Fol- 
lowing Ref.[46jone can introduce a unitary matrix 9 which 
satisfies the equation: 

wr6» t = 0tV = diag{VTU, 

where all quantities are energy dependent. The solution 
is 6 = diag{l/v / 7^}wT. The following properties of 9 can 
be checked: 99^ = 5 VV <, 9^9 = S kk ', thus 9 is indeed the 
unitary matrix. It diagonalizes rV: 



0( T t T )0t = diag{T y }. 



(6) 



Matrix 9 performs a unitary transform of out-states in R, 
so that the linear combination v of in-coming states 
in L scatters into the linear combination v of the out- 
states in R, $° out (r, E)\ z ^ +00 = £ fc 6 v k (E) $ £ out (r, E), 
with the transmission amplitude ^/T V (E), namely: 

*„(r, £)U +oc = y/TjE) 4>° out (r, E). 

One can sho w 28 ' 31 that for the Bloch states at the same 
energy (Ek = E = E k ') in the ideal leads the following 
relations hold for the current matrix elements: 



dS 



%(r,E)id z %f(r,E) 



f^<W, (7) 



where the Bloch waves are either left (v K oc (dEk/ '9k) 2 > 

0) or right-travelling (v K < 0), the operator d z is defined 

as / dz 9 — f(d z g) — (d z f)g, and the integral goes over 
infinite plane S (cross-section of the lead) which is per- 
pendicular to the current direction z. 

In case of the perturbed system the orthogonality re- 
lation holds for current matrix elements in the basis of 
eigenchannels. Using Eqs. ([3]), ([5]) and ([7]) we can com- 
pute it in the asymptotic region of the right (R) lead: 47 



dS 



<S>„(r,E)id z K(r,E) 



z— *+oo 



ME) 
2tt 



(8) 



We note, that Eq. ([5]) holds for any position z of the 
plane S. Because the wave functions of channels are so- 
lutions of the Schrodinger equation with a real potential 
corresponding to the same energy, the flux through arbi- 
trary plane S is conserved (a proof is similar to that of 
Appendix A of Ref. \3lh . 



IV. CONDUCTION EIGENCHANNELS WITHIN 
THE KKR METHOD 



&v(r,E) 



(4) 



A. Evaluation of conductance 



The unitary transform u> is defined such way that the 
transmission matrix T — tt^ is diagonal in the basis v. 



w(Trt)o;t = diag{T I/ }, 



(5) 



and the conductance reads as g = go^2, v T v (Ep), where 
the T„'s are transmission probabilities of eigenchannels. 



To calculate the ballistic conductance, we employ the 
Kubo linear response theory as formulated by Baranger 
and Stoned 

9 = 9o / dS dS' 

x G + (r,r',E F ) d z d z > G-(r',r,E F ), (9) 



where G~ and G + are retarded and advanced Green's 
functions, respectively, and go = 2e 2 /h. The integration 
is performed over left (Sl) and right (Sr) planes which 
connect the leads with the scattering region (Fig. 1). 
The current flows in z direction. The implementation 
of Eq. |9]) within the KKR method, related conver- 
gence tests and further details were discussed in recent 
publication^ 

In this paper we present a further extension of Ref. l3ll . 
a method for the evaluation of conduction eigenchannels. 
For that, we will follow closely the analysis of Refs. 28,31 
where the equivalence of the Kubo and Landauer ap- 
proaches to conductance problem was shown. We pro- 
ceed in four steps: (i) we remind the KKR representa- 
tion of the Bloch functions; (ii) we build up the asymp- 
totic expansion of the Green's function in terms of un- 
perturbed Bloch states of the leads; (iii) we construct 
the transmission matrix T = tt^ in fc-space; (iv) finally, 
we find an equivalent representation of the transmission 
matrix within the local KKR basis. Further solution of 
the eigenvalue problem leads us to conduction channels. 
We mention here one aspect of the problem: for realistic 
calculations the planes Sl and Sr (see Fig. 2 for details) 
are usually placed in finite distance from the atomic con- 
striction. Nevertheless, we first focus on the asymptotic 
limit before we discuss the realistic situation which is 
considered in Sec. IV. E and Appendix A. 



B. Atomic orbitals and Bloch functions 

Let r be arbitrary point in the asymptotic region of 
the lead (Fig. 2). In the KKR method the local, energy 
dependent basis of atomic orbitals is defined at each unit 
cell n: 



supercell V L 



<%(t,E) = L (r-R„,£) 

= R L (r-TL n ,E)<d(r~R n ), 



(10) 



where Rl's are real regular solutions^ of the Schrodinger 
equation for the potential V n (r) at cell n, and 0-function 
is 1 inside cell n and is outside it. The unperturbed 
Bloch function is given by expansion over atomic orbitals 
at all sites n in the Born-von Karman supercell: 

n(r,E) =Y,C k ME)Mr-'Rn,E) (11) 

iiL 

with C k , nL (E) = e+ lkR "C£ L {E). Here the common 
notation k = (k, A) for the Bloch vector k of the 
1st Brillouin zone (BZ) and the band index A is used. 
The C kL {E) are solutions of the KKR band structure 
equations^ with energy E. 

Considering a waveguide geometry, we will assume the 
Bloch functions to be normalized per cross section of 
the Born-von Karman supercell (see Fig. 2) with open 
boundary conditions along z-axis. Thus, Ck,nh{E) oc 
1/ yjN x N y with N x N y being number of atoms per cross- 
section, and orthogonality condition for Bloch waves 
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FIG. 2: Sketch of the system under consideration. In a for- 
mal theory, left (£>£) and right (S R ) planes are placed within 
the leads far away from a scattering region (asymptotic limit). 
For all points r (z < z£) and r' (z! > z^) within the Born-von 
Karman supercells (cubes) Vl and Vr asymptotic properties 
are achieved. When implementation of the method is con- 
sidered, the conductance is calculated between the planes Sl 
and Sr positioned somewhere in the scattering region. 



takes a form 31 : 



J d 3 rn(r,E)^° k r(r,E') = 



V 



-S 



fcfc' 



(N X Ny) 

v 

= \v k \6{E-E')5 ss ,5 HkV 

Here V — L x L y L z is a volume of the supercell, Aq = 
L x L y /N x N z is an area of the xy-xmit cell in the electrode, 
s,s' = ±1 are signs of k z and k' z , relation L z Sk z k' = 
2-7T S(k z —k' z ) has been used, and velocity v k along current 
flow is defined as v k — ZnA^dEk/dtyz. 

Since the Bloch waves form a complete set, a back 
transform of Eq. (fTT|) exists: 

k' 

where k' sum runs over all k'-points in the 1st BZ and 
over all bands A'. The f symbol means Hermitian conju- 
gate. The expression for B\ L k can be obtained^ from 
known matrix C k ,nL- One can prove further, that C kiTl L 
and -B„ L k obey the following orthogonality relations: 



E 

ii L 



C k , nL (E) Bl L k ,(E) = <W, 

Y, B nL, k (E)C k , n/L/ (E) = S nn ,5 LL/ , (13) 
k 

where in the second equation a sum over k is restricted 
to states with E k = E. 

C. Asymptotic expansion of the Green's function 

Starting from the site angular momentum represen- 
tation ([T]) of the retarded Green's function within the 
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KKR method and using Eqs. (jlip and (fT2|) we obtain 
the asymptotic expansion for G + (r,r',E) over the un- 
perturbed Bloch waves ! 28 ' 31 



z,z' -H>=poo 



(14) 



kk' 



with r e V L , r' e Vr (see Fig. 2), and 

A kk >{E)= (15) 
E E E^."^^^^)^,,,,,^) 

nSVL n'GVR LL' 

or in a matrix form: A = BGB X . Formally, the fc-sums 
in Eq. (TT3]) are performed over all k-states in the 1st BZ 
and over all bands A. However, since the Green's func- 
tion for r ^ r' is a solution of the Schrodinger equation 
without a source term, only states fc, fc' at the isoener- 
getic surface of energy E contribute to the sum in the 
asymptotic expansion! 28 ' 31 Therefore, 

A kk '(E) = -^5(E-E k )5(E-E k ,)A kk ,(E), (16) 

^ Z 

where S(E - E k ) = (fi /S A ) Wk\ 5{E - E k ). Here N z 
is number of atom sites in Born-von Karman supercell 
along z axis, Slo = (2ir) 3 /Vo is volume of the 1st BZ, 
S\ is the area of the isoenergetic surface corresponding 
to band A, v k = dE k /d]s.. For the discrete fe-points, 
the function (l/N z )S(E — E k ) equals 1 if E = E k , and 
is otherwise. In addition, boundary conditions for the 
Green's function 2 ^, constrain matrix elements A kk > (E) to 
be non-zeros only if k and fc'-states are right-travelling 
waves (with positive velocity along z-axis) that corre- 
sponds to the in-states k in the left lead and out-states 
k' in the right one. 



D. Transmission matrix: asymptotic limit 

We proceed further, and use the asymptotic represen- 
tation (TT4"|) of the Green's function to evaluate conduc- 
tance according to Eq. ([9]). Assuming the integration 
planes to be placed within the leads infinitely far from 
the scattering region, we obtain: 



9 



g Tr {k) [VlA(E f )V k A1{E f )] , 



(17) 



where the diagonal operators of velocities Vl and Va (re- 
lated to the left and right planes) acting in the fc-space 
were introduced: [Vl(r)] fefe , = (v k /2ir) 5 kk >. Formally, 
the trace (Tr) in Eq. (fTTj) goes over all fc-states and the 
Fermi surface is taken into account by means of Eq. (|16p 
where E = Ep. 

The velocity operators can be decomposed into sum 
of two operators related to the Bloch states with pos- 
itive and negative velocities along z: V = V + + V~ , 



where V + is nonzero for right-travelling waves only, while 
V~ is nonzero for left-travelling ones. In the asymptotic 
limit, only in-coming and out-going fc-states with posi- 
tive velocities contribute to the sums in Eq. p"?)) . Using 
the relation between expansion coefficients A k y and the 
transmission amplitudes r kk > derived in Refs.TSlll 



ifcfc 



(E) 



-2m 



Tkk'(E) 



we obtain: 



g Ti (h) [V+ A(Ep) V+ A^ (E F )] 
ffoTr (fc) [T(E F )], 



(18) 



(19) 



where a representation of T = tt^ in fc-space is given by 

T(E) = {V+) l ' 2 A{E) V+A^E) {V+) 1 ' 2 (20) 

with a positive definite operator under square-roots. 

The fc-representation is formal but not suitable for 
implementation. To solve the eigenvalue problem for 
T = tt^ the mapping on the site-angular momentum 
(n, L)-space of the KKR method should be presented. 
Such mapping is realized through the expansion (fTTj) of 
the Bloch functions over atomic orbitals, so that velocity 
operators in fc-space take a form: 

[ V H*)]kk' = E 22 C k,nL [ D HR)]TL' C n'L',k" 
nn'eSl LL' 

where site-diagonal operators -Dr and £>l defined on 
atomic orbitals are the KKR-analogue of velocities^ 



[D H r)(E)] 

±5nn' 



LL' 



(21) 



dS 



Rl(r,E)id z Rl,(r,E) 



here integral is restricted to the cross-section of the unit 
cell around site n. Now we can evaluate conductance 
according to Eq. (fTT)) . Taking into account that B^C — 
C^B = 8 nn t5 L Lt [Eq.(rj3])], we obtain: 



a 



= <?oTr ( „ iL) (D L GD R G^) 



(22) 



here (and further) all matrices are assumed to be taken at 
the Fermi energy, G = {G"^'} stands for matrix notation 
of the structural Green's function introduced in Eq. |T]), 
and trace (Tr) involves sites and orbitals related to the 
atomic plane 5° (Fig. 2). 

The operators D\ J1 Dr are anti-symmetric and hermi- 
tian, thus their spectrum consists of pairs of positive and 
negative eigenvalues: ±<5°. Let U be the unitary trans- 
form of matrix D (cither or Z?r) to a diagonal form: 



D° = U^DU = diag{±<5°}. 



(23) 



Decomposition of operator D into two terms, D = D + + 
D~ , related to the right- and left-travelling waves is nat- 
urally given in the basis of eigenvectors: 



D ± = UD o{±)u\ 



(24) 



6 



with 



A+ 




D o(- 





-A" 



where A + , A are positive (non- negative) diagonal ma- 
trices. The analogue of Eq. (fT9"|) reads as 



9\ 



{D+GD+ Gt) 



(25) 



Now we are ready to build the (n, L)-representation 
for transmission matrix T — tt< , For that, one should 
extract a square-root from the positive definite operator 
Vl = CD^C' defined on fe-space with help of (n,L)- 
space. Namely, because the is positive definite, it 
can be represented in the following form: 52 = SIl^l- 
Here an operator f^L maps the fe-space on the (n, L)- 
space, J7l : {fe} — * {n, L}. The solution for Qj, is 

n h = C(D+) 1 / 2 £, 

where £ is an arbitrary unitary matrix (££^ = 1) and 

{D+) 1 ' 2 = U L [Dl i+) ]^ 2 Ul with [£ L (+) ] 1/2 being square- 
root of the positive definite diagonal matrix [Eq.([24])]. 

To find T = tt\ we start from Eq. ([19]) and proceed 
as follows: 



5oTr (n,L) (nlAv+An^ 

~{D+y> 2 C\BGB^){CD+&) 



= ffoTf(n,i) 



{BG^B^)C{Dl 



a/2 



where Eq. (fl~5)) was used. Because of B^C — 5 nn t5LL> 
[Eq. ([13])]. we obtain: 



9 = 90^(n,L)(T), 



(26) 



where, in the asymptotic limit, the transmission matrix 
T = tt^ in (n, L)-representation is given by 



T={Dt)^GD+G\D+Y 



/2 



(27) 



The trace in Eq. ([26)1 goes over all sites n and orbitals L of 
the atomic plane 5° in the left lead (Fig. 2). Equivalent 
formula can be written for the right lead. To conclude, 
one can prove that spectrum of the obtained matrix coin- 
cides with spectrum of transmission matrix T [Eq. ([2H)l ] 
defined in fe-space (see Appendix B for details). There- 
fore, solution of the eigenvalue problem for T gives us 
required transmission probabilities. 



E. General case: arbitrary positions of planes 

In practical calculations of conductance with the use 
of the Kubo formula, the left (L) and right (R) planes are 
positioned somewhere in the leads (Fig. 2). Expression 
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FIG. 3: Top (a): A model step-like potential for the free- 
electron gas within the L/S/R structure having two dimen- 
sional periodicity. The electrons moving with energy E — Ep 
from the leads (L and R) are scattered at the potential 
step U. The fey momentum is conserved, while the momen- 
tum k z along z-axis is ko — \2mE — fc|] ly/2 in the leads, and 
fci = [2m(E - U) - fcjj] 1/2 in the S layer. Bottom (b): A mul- 
tiple reflection of the in-coming from the R lead wave e~ lk ° z 
within the spacing layer S, where it is a linear combination 
4>{z) = x+(z)+X-(z) = ae~ ihiz +/3e +lhiz of two functions. 
If z L is position of one of the planes taken for the conductance 
evaluation, transmission T(k«) of channel feii is proportional 

to the current j = (j)(z) id z <j)*(z) = 2k = 2ki \a\ 2 - 2ki \/3\ 2 = 
j+ — j- , which is sum of two terms (positive and negative) 
due to contributions from two functions, x+( z ) ana - X-( z )- 
The value p = f3/a = (fci — feo)/(fei + k )e +2lkia has a mean- 
ing of the reflection coefficient of the S/R or S/L interface. 
However, if plane is chosen in the asymptotic region of the 
lead, at point z®, transmission T(fcii) ~ 2ko contains only one 
positive contribution. 



is valid in general case and result is exactly the same 
as in Ref. Operator in Eq. ([2"2"l) is sum of two 
terms: 



D L . Therefore, we can write down 



g = g Tr(D+GD R Gl 

= <7 + + <T, 



go 



Tr (D^GD K G"<) 



(28) 



where g + and g denote two contributions. In a formal 
theory, when the atomic plane 5l is placed in the asymp- 
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totic region of the left lead far from the atomic constric- 
tion the second term g~ in Eq. (f2"8|) is equal to zero. In 
practice, the real space summation of current contribu- 
tions includes only a finite number of sites at both atomic 
planes, because the current flow along z direction is lo- 
calized in the vicinity of the contact. Due to numerical 
effort we are forced to take integration planes closer to the 
constriction in order to obtain convergent value for the 
conductance with respect to number of atoms included 
in summation. In addition, even better convergence for 
matrix elements is required to solve the eigenvalue prob- 
lem. A compromise can be usually achieved but positions 
of the atomic planes Sl and Sr do not meet the asymp- 
totic limit criterion. However, since the electron current 
through the structure is conserved, any position of the 
planes is suitable for the calculation of conductance. If 
Sl is placed somewhere in the scattering region we have 
to sum up all multiple scattering contributions. We show 
in Appendix A, that all multiple scattering contributions 
in direction of the current cause g + , whereas all scatter- 
ing contributions in opposite direction give rise to g~ . 
Thus, the first term, g + , in Eq. (I28|) is always positive, 
while the second one, g~ , is always negative. To make 
this statement clear, an illustration of scattering events 
is shown in Fig. 3 assuming a simple free-electron model. 
In the region of the lead where the potential is a small 
perturbation with respect to the bulk potential the con- 
tribution to the conductance due to g~ is one order of 
magnitude smaller than g + . 

To find transmission probabilities of eigenchannels, one 
has to apply the procedure introduced in the previous 
section independently for both terms, g + and g~ . We re- 
fer to Appendix A for a mathematical justification. Ex- 
pression for conductance takes a form: 

9 = .9oTr ( „, L) (T+) + g Tr („,/,) (T~), (29) 

with 

T± =±{±D±) 1/2 GD K G\±D±f' 2 . (30) 

We show in Appendix A that all eigenvalues of T + are 
either positive or zeros whereas all eigenvalues of are 
negative or zeros. To identify transmission probabilities 
T n of channels the spectra of operators T + and T" have 
to be arranged in a proper way. Then transmission of the 
rt-th channel is given by T n = t+ — t~ , where irj^ are 
positive and negative eigenvalues of the operators T + 
and T~, respectively. The T„ does not depend on the 
positions of the left (Sl) and right (Sr.) planes, while 
are z-dependent. In the asymptotic limit t~\ Zl ^^qc — ► 
and T n = t£\ Zl ^-oo, so that the Landauer picture is 
restored. 

In general case the direct way to find the pairs of eigen- 
values is not evident without a back transform to the k- 
space. However, from the point of view of applications 
to the extremely small symmetric atomic contacts, as the 
ones we are studying in this work, the problem is easy to 
handle. Since the number of contributing eigenmodes is 



limited, the pairs of eigenvalues can be found by symme- 
try analysis of the eigenvectors of T + and T~ . Namely, 
using the symmetry properties of the structural Green's 
function GYl> and current matrix elements D^, one can 
show that channel's transmission T n bounded between 
and 1 is defined by eigenvalues t+ and t~ which belong 
to the same irreducible representation of the symmetry 
point group. 



V. APPLICATIONS OF THE METHOD 

In recent paper o 29 ' 30 we have verified the method de- 
scribed here by studying systematic changes in the con- 
ductance of metallic constrictions in the presence of de- 
fect atoms. Illustrative examples presented below focus 
on single-atom contacts made of pure metals such as Cu, 
Ni, Co and Pd. 

Copper serves mainly for test purposes. It is a repre- 
sentative of the noble metals and has electronic prop- 
erties similar to Ag and Au, for which a lot of ex- 
perimental results^ 3 ^ 3 -! 5 ^ 5 ^ 6 - as well as DFT based 
calculation s 21 ! 24 are available. In particular, a large num- 
ber of experiments for alkali metals (Li, Na, K) 32 i 56 i 57 
and noble metals (Au, Ag, Cu); 9 - ' 32 ' 54 ' 55 ' 56 employing dif- 
ferent techniques under room an liquid-He temperatures, 
show that conductance histograms have a dominant peak 
very close to one conductance quantum Go = 2e 2 /h, and 
smaller peaks close to integer values. 

However, for transition metal contacts (examples 
of which are Ni, Co and Pd) the situation differs 
significantly . 32 ' 55 Only one broad maximum centered 
somewhere between I Go and 3 Go is usually observed 
in conductance histograms. 4 ^ That is a signature of 
the nontrivial decomposition of conductance consisting 
of more than one perfectly transmitting channel ) 15 ' 16 
since for transition metal atoms d states of different 
symmetry are available at the Fermi level. The ques- 
tion on half-integer conductance quantization has been 
addressed i 33 ' 58 However, recent experiments^ do not 
confirm this hypothesis thus pointing to the conclusion 
that the electron transport through ferromagnetic con- 
tacts can never be fully spin-polarized. Another issue is 
a large magnetoresistance (MR) effec t 40 ' 41 ' 42 observed for 
metallic point contacts made of different magnetic mate- 
rials. This field is known to be full of controversy. There 
is still a continuing discussion on whether or not the enor- 
mous MR values found experimentally could be of elec- 
tronic origin^i 3 ^ 6 ^!^! 6 ^ 3 .!^ 4 . or the effect is just due to 
atomic rearrangements in the neck region of a contact as 
a response to the applied magnetic field i 65 ' 66 Extensive 
discussion on this topic can be found in the recent re- 
view paper by Marrows^i We will comment further on 
the above issues. 
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FIG. 4: (Color online) Energy- dependent transmission and 
its decomposition to the conduction eigenchannels for the Cu 
single atom contact shown in Fig. 1. 



A. Computational details 

An atomic configuration of a constriction used in our 
calculations is presented in Fig. 1. The single atom con- 
tact was modeled by a small cluster attached to the semi- 
infinite fee (001) leads. The cluster consists of two pyra- 
mids joined via the vertex atoms separated by a distance 
. Here oo is the experimental lattice constant of 
fee metals: 6.83 a.u. for Cu, 6.66 a.u. for Ni, 6.70 a.u. for 
Co, and 7.35 a.u. for Pd. The metals under considera- 
tion do not have a tendency to form chains.— An atomic 
bridge is most likely to be broken just after the single 
atom limit is achieved. Thus, a configuration shown in 
Fig. 1 resembles one of limiting configurations of point 
contacts which could appear in the MCBJ experiments. 

Our calculations are based on DFT within the local 
density approximation. The parametrization of Vosko, 
Wilk, and Nusair— for the exchange and correlation en- 
ergy was used. The potentials were assumed to be spher- 
ically symmetric around each atom (atomic sphere ap- 
proximation, ASA). However, the full charge density, 
rather than its spherically symmetric part, was taken 
into account. To achieve well converged results the an- 
gular momentum cut-off for the wave functions and the 
Green's function was chosen to be l max — 3 that imposed 
a natural cut-off 2/ roax = 6 for the charge density expan- 
sion. In case of heavy element Pd the scalar relativistic 
approximation^ was employed. For the conductance cal- 
culation the surface Green's function was computed us- 
ing a small imaginary part ImE = 0.04 mRy and about 
250.000 k-points were taken in the 2D Brillouin zone. In- 
stead of integration over planes, current matrix elements 
(|2"Tj) were averaged over atomic layers as described in 
detail in Ref. [3l|. A typical error in the calculation of 
conductance was ~ 5%. 



B. Symmetry analysis of eigenchannels 

To understand the relation between the electronic 
structure and the transport properties of atomic con- 
tacts we consider the energy-dependent total transmis- 
sion, T(E), and its decomposition to the conduction 
eigenchannels, Ti(E). Results are shown in Figs. 4, 6, 8, 9 
for the case of Cu, Ni, Co and Pd point contacts, re- 
spectively. The investigated structure (Fig. 1) has a C± v 
symmetry. Further we denote individual channels by the 
indices of irreducible representations of this group using 
notations of Ref. [7ll. common in band theory. In addition, 
each channel can be classified according to the angular 
momentum contributions when the channel wave func- 
tion is projected on the contact atom of the constriction. 
This is very helpful since the channel transmission can 
be related to the states of the contact atomJ^ For exam- 
ple, the identity representation Ai of the C± v group is 
compatible with the s, p z and d z i orbitals (here the z is 
perpendicular to the surface and passes through the con- 
tact atom), while the two-dimensional representation A5 
is compatible with thepz, p y , d xz , d yz orbitals. The basis 
functions of A2 and A2' are d x 2_ y 2 and d xy harmonics, 
respectively. 



C. Cu contacts 

The energy-dependent transmission of Cu atomic con- 
tact (shown in Fig. 1) is presented in Fig. 4 together with 
the eigenchannel decomposition. At the Fermi energy the 
calculated conductance value is G — 1.01 Go- It mainly 
consists of one open channel of Ai symmetry which arises 
locally from s, p z and d z 2 orbitals when the wave func- 
tion is projected on the contact atom. This result is in a 
good agreement with a lot of experiment s 32 ' 56 mentioned 
previously as well as with other calculations involving dif- 
ferent approaches! 21 ' 24 The additional twofold degenerate 
channel has A5 symmetry. Transmission of this channel 
increases at energies above the Fermi level (Ep) together 
with an increase of the p x ,p y contribution to the local 
density of states (LDOS) at the contact atom. However, 
at E < —1.5 eV below the Ep the A5 channel is built 
mainly from the d xz , d yz orbitals of the Cu atom. 

We would like to point out that in case of noble metals 
the conductance of single-atom contact is not necessarily 
restricted to one channel. An example of a configura- 
tion which has more channels (but still has only one Cu 
atom at the central position) is presented on the top of 
Fig. 5. Because of the larger opening angle for incoming 
waves as compared with the preceding case, conductance 
of such system is G = 2.57 Go with major contribution 
from four channels (Fig. 5 and its caption). The value 
2.57 Go correlates with a position of the third peak in the 
conductance histogram of Cu, which is shifted from 3 Go 
to smaller valuesj 32 ' 56 As it is seen from the presented ex- 
ample, conductance quantization does not occur for the 
metallic atomic-sized contacts. In general, even for no- 
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majority spin 
minority spin 



channel 1: A 



S-2 plane 

channel 4: A. 




channel 2: A, 



channel 3: A, 





FIG. 5: Wave functions |^i(r)| 2 (probability densities) of the 
four dominating eigenchannels for the pyramidal Cu contact 
shown on the top. Wave functions resolved to atoms are visu- 
alized 2 atomic planes below the surface plane (S-2). Colors 
from white to black correspond to consequently decreasing 
positive values. Transmission probabilities of channels are: 
Ti = 0.90 (Ai), T 2 = T 3 = 0.71 (A 5 ), T 4 = 0.08 (Ai), which 
are summed up to conductance G = 2.57Go- Further details 
are given in the text. 



ble metals, conduction channels are only partially open 32 
in contrary to the case of quantum point contacts real- 
ized in the two-dimensional electron gas where a clear 
conductance quantization was observed.— 

For illustration, we present in Fig. 5 probability ampli- 
tudes of the eigenchannles |\&i(r)| 2 in real space resolved 
with respect to atoms of the 2nd plane (S-2) below the 
surface. We see that the wave functions of the 1st and 
the 4th channel with the highest symmetry (Ai) obey all 
eight symmetry transformations of the C± v group, while 
two different wave functions of the double degenerate 
channel (A5) are transformed to each other after some 
symmetry operations. 



D. Transition metal contacts 

We turn to transition metals, and consider the ferro- 
magnetic Ni assuming a uniform magnetization of the 
sample. Transmission T(E) split per spin of a Ni contact 
is shown in Fig. 6. A shift about ~ 0.8 eV along the 
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FIG. 6: (Color online) Spin-dependent transmission as a func- 
tion of energy for the Ni single atom contact shown in Fig. 1. 



energy axis between transmission curves is seen that is 
in agreement with exchange-splitting of the Ni d states. 
Similar computational results regarding transmission of 
Ni constrictions were reported by Solanki et al., 36 Rocha 
et alw^ and Smogunov et al.£2- Exchange splitting esti- 
mated from their works varies from 0.8 eV till 1.0 eV, 
but fine details differ because of different atomic config- 
urations and employed methodologies. In this regard, 
exchange splitting about 2.0 eV observed in transport 
calculations of Jacob et al. 37 in case of Ni contact seems 
to be overestimated. 

The shift in energy due to different spins is observed as 
well for the transmissions of individual channels (Fig. 7). 
We see from Fig. 7 and Table I that at the Fermi en- 
ergy the spin-up (majority) conductance of Ni contact 
is mainly determined by one open Ai channel (similar 
to the case of Cu), while three partially open channels, 
of Ai and A5 symmetry, contribute to the spin-down 
(minority) conductance. The minority A5 channel arises 
locally from d xz and d yz states, rather than from p x ,p y 
states whose contribution to the spin-down LDOS at the 
Fermi energy is much smaller (Fig. 7). The calculated 
conductance, G — 1.20 Go, correlates with a position of 
the wide peak in the conductance histogram of Ni cen- 
tered between 1 G and 2 G &&- 

Within the energy range shown in Figs. 6 and 7 (±2 eV 
around Ep), we count six eigenmodes of different sym- 
metry for both spins. At energies well above the Fermi 
level (E > 1.0 eV) the spin-splitting of sp Ni states is 
lost, and the picture is similar to what we have seen 
for Cu. Three channels are present: one open Ai (sp z - 
like) channel with transmission around 0.9 and a partially 
open double degenerate A5 (p Xl p y ) channel whose trans- 
mission increases monotonically as a function of energy. 
However, below the Fermi energy all eigenmodes Ti(E) 
display a complicated behavior (Fig. 7, upper plots) that 
reflects a complex structure of the LDOS projected on 
orbitals of the contact atom (Fig. 7, bottom plots). Be- 
low Ep the existing s, p z and d z 2 states are strongly 
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FIG. 7: (Color online) Spin- and energy-dependent transmission decomposed to conduction eigenchannels for the Ni atomic- 
sized constriction in comparison with the symmetry projected local density of states at the contact atom (i.e. apex atoms in 
Fig. 1). 



TABLE I: Transmission probabilities of eigenchannels at the Fermi energy of Ni and Co atomic contacts shown in Fig. 1 for two 
different (P and AP) orientations of magnetizations in the leads. Only transmissions of the dominant channels are presented. 
Magnetoresistance ratio defined as MR = (Gp — Gap) I Gap x 100% is given in the last line. 







Ni 






Co 




Channel 


P 




AP 


P 




AP 




CD-spin 


(T)-spin 


(T) or (I) spin 


(|)-spin 


(t)-spin 


(T) or (I) spin 


Ti (Ai) 


0.68 


0.84 


0.82 


0.36 


0.89 


0.58 


T 2 = T 3 (A 5 ) 


0.35 


0.06 


0.31 


0.14 


0.07 


0.09 


T 4 (A 2 ) 








0.17 






Transmission 


1.44 


0.96 


1.45 


0.83 


1.03 


0.76 


MR ratio 




-17% 






+23% 
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hybridized giving rise to two channels of Ai symmetry 
at energies about E m —1.0 eV and E » — 1.5 eV for mi- 
nority and majority spins, respectively. A clear correla- 
tion between the symmetry projected LDOS and Ti(E) is 
seen for the pure d-channels, A2 (d x 2_ y 2) and A2' (d xy ). 
For example, the minority spin d x 2_ y s resonance centered 
around E = —0.7 eV and majority spin d x 2_ y 2 resonance 
at E = —1.3 eV reflect themselves as peaks in the trans- 
mission of the minority and majority A2 channels. The 
same is valid for the d xy states at E = —1.1 eV (spin 
down) and E = —1.7 eV (spin up) which cause the in- 
crease of transmission of the A2' channel at the same 
energies. However, even at the resonances, the A2 and 
A2' channels are only partially open because the d x 2_ y 2 
and d xy orbitals are spread perpendicular to wire (z-) 
axis that prevents effective coupling with the neighbor- 
ing atoms. 

Our results for Co constriction are presented in Fig. 8. 
As compared with Ni, the shift between the spin-up and 
the spin-down T(E) curves becomes larger (~ 1.7 eV), 
because of the stronger exchange field of cobalt. At the 
Fermi energy majority spin conductance is still domi- 
nated by one highly transmitted Ai channel (Table I), 
while for the minority spin the d x 2 _ y 2 resonance is pinned 
to the Fermi level and results in the additional (as com- 
pared with Ni) channel of A 2 symmetry. Thus four chan- 
nels with moderate transmission probabilities contribute 
to the minority spin conductance (Table I). 

Fig. 9 shows results for Pd. According to recent the- 
oretical predictions^ 3 - monatomic Pd wires might exhibit 
magnetic properties. However, in this work we considered 
nonmagnetic solution, since a coordination number even 
for the contact atom was already big enough (Fig. 1) to 
suppress magnetism. 73 Pd is isovalent to Ni. The Fermi 
level crosses the partially filled d band. Therefore, the 
eigenchannel decomposition resembles the minority spin 
channels of Ni. However, due to a larger occupation num- 
ber, transmission curves are shifted ^0.5 eV downwards 
in energy as compared with spin-down Ni modes. The 
conductance G — 1.41 Go is a sum of three channels. This 
value is in agreement with the conductance histogram of 
Pd2i which shows a broad maximum around ~ 1.7 Go- 

We turn back to Ni and Co contacts and consider a sit- 
uation when a relative orientation of magnetizations in 
the leads is antiparallel (AP), so that an abrupt atomic- 
scale domain wall is formed as shown in Table II. We see 
from Table I that, both for Ni and Co, the AP conduc- 
tance reflects the structure of the minority spin channels 
and consists of a Ai and a A5 channel. For the atomic 
configuration shown in Fig. 1, we obtained "optimistic" 
MR values: —17% in case of Ni, and +23% in case of Co, 
which arc quite small in accordance with our previous 
studyi^ 5 - However, the precise MR values as well as trans- 
mission curves for Co, differ from the results reported in 
our earlier work because of the different geometrical con- 
figurations of atomic contacts. The reason is that the 
transmission of d-like channels is quite sensitive to the 
exact geometry* 3 ^ We mention here, that a more accu- 
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FIG. 8: (Color online) Spin- and energy-dependent transmis- 
sion decomposed to conduction eigenchannels for the Co sin- 
gle atom contact shown in Fig. 1. 



rate full-potential approach and an improved description 
of the electron correlations for localized (i-electrons can 
somewhat affect presented results. That is also valid for 
the effects of atomic relaxations which were neglected. In 
particular, the exact values for the transmission proba- 
bilities and MR at the Fermi level reported in this study 
for different systems could be slightly changed. How- 
ever, more precise calculations obviously will not affect 
the physical results of the present work. 

Evident conclusions follow from presented examples. 
First, in contrast to earlier studies , 33 ! 58 the ferromagnetic 
Ni and Co contacts do not show any tendency to close 
one spin channel. On the contrary, both spin channels 
contribute to the conductance that gives only moderate 
magnetoresistance values. Independent of the geometry 
of the atomic contact the minority spin channel will in- 
clude a sum of fractional contributions from many modes 



12 



UJ 



c 
o 

in 
</> 

E 
(A 
c 
re 




-Total transmission 
-channel 1: A 

-channels 2,3: A 5 

-channel 4: A 2 

channel 5: A 

-channel 6: A, 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 
£-£ F [eV] 

FIG. 9: (Color online) Eigenchannel decomposition of the 
transmission for the Pd contact shown in Fig. 1. 



TABLE II: Spin magnetic moments (in /is) at atoms form- 
ing Ni and Co contacts shown in Fig. 1 for the parallel (P) 
and the antiparallel (AP) orientation of magnetizations in 
the leads. Bulk magnetic moments are 0.62 /i_g for Ni, and 
1.62 ub for Co. 



Ni 



Co 



Atom 


P 


AP 


P 


AP 


Surface 


0.66 


0.66 


1.78 


1.78 


Contact— 1 


0.70 


0.69 


1.83 


1.83 


1st contact 


0.68 


0.54 


1.76 


1.64 


2nd contact 


0.68 


-0.54 


1.76 


-1.64 


Contact+1 


0.70 


-0.69 


1.83 


-1.83 


Surface 


0.66 


-0.66 


1.78 


-1.78 



because the d states are always present at the Fermi level. 
That agrees with later experiments by Untiedt et al.^ 
where the absence of conductance quantization for ferro- 
magnetic Fe, Co and Ni contacts was clearly confirmed. 

Second, an abrupt, atomic-scale domain wall pinned 
to the constriction does not show an impressive MR ef- 
fect. For a fixed atomic configuration, the P and AP 
conductances are of the same order. Most likely, that 
more sophisticated calculations, involving relaxation ef- 
fects and noncollinear magnetic moments in the domain 
wall, will not be able to change this statement / 5 ' 76 A re- 
cent research 7 - 7 - towards transport in nanocontacts with 
noncollinear moments shows that energetically preferable 
noncollinear magnetic order results in a larger domain 
wall width as compared to the abrupt, collinear wall con- 
sidered in the present paper. That leads to weakened 
scattering of electrons and a further reduction of the MR 
values. 

Coming to the experimental situation on ballistic mag- 
netoresistance (BMR) effect in ferromagnetic contacts, 
we point out that large MR value s 41 ' 42 were usually mea- 



sured for much thicker constrictions (as compared with 
atomic-sized contacts) with resistance in the range of 
hundreds of Ohms. It is believed ) 65 ' 66 that such exper- 
iments suffer from many unavoidable artifacts induced 
by magnetomechanical effects that mimics the real MR 
signal which would come from the spin-polarized trans- 
port alone. However, recent studies by Sullivan et 
and Chopra et alM on Ni and Co atomic-sized contacts 
report BMR values in the range of 200 4- 2000% with 
discussion on the electronic origin of the effect due to 
domain wall scattering. Nevertheless attempts to min- 
imize magnetostrictive effects were undertaken, we just 
can repeat^ 5 - that a natural explanation of thes o 43 ' 44 and 
similar experiment o 34 ' 40 is that, due to magnetization 
reversal processes, unstable in time atomic constriction 
changes its contact area when magnetic field is applied. 
Characteristic steps and jumps in the measured field- 
dependent conductance (Fig. 4 of Ref. I43T ) or resistance 
(Fig. 3a of Ref. 0, Fig. 3 of Ref. Hi) are distinct evi- 
dence of atomic reconstructions and fractional changes of 
the contact cross section. For example, just eliminating 
one contact atom from the configuration shown in Fig. 1 
changes conductance of a Ni constriction from 1.2 Go 
(chain of two atoms, see Table I) up to ~ 2.8 Go (one 
contact atom only, see Ref. [13), thus producing ~ 130% 
MR. Further increase of a contact area can give arbitrary 
high MR values, that supports hypothesis on mechanical 
nature of the effect. 



VI. CONCLUSIONS 

To summarize, we have presented a formalism for 
the evaluation of conduction eigenchanncls of metallic 
atomic-sized contacts from first-principles. We have com- 
bined the ab initio KKR Green's function approach with 
the Kubo linear response theory. Starting from the scat- 
tering wave formulation of the conductance problem, we 
have built a special representation of the transmission 
matrix in terms of local, energy and angular momen- 
tum dependent basis inherent to the KKR method. We 
have proven that solutions of the eigenvalue problem for 
the obtained matrix are identical to conduction eigen- 
channels introduced by Landauer and Buttiker. Appli- 
cations of the method have been presented by studying 
ballistic electron transport through Cu, Pd, Ni and Co 
single-atom contacts. The symmetry analysis of eigen- 
channels and its connection to the orbital classification 
known from the tight-binding approach were discussed 
in detail. Experiments on the electron transport through 
magnetic contacts were commented. 
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APPENDIX A 

We consider a practical implementation of the Kubo 
formula ([9]): positions z = zl and z' = zr, of the Sl and 
Sr planes are chosen within the semi-infinite electrodes, 
where Ai-matrix describing a scattering region is suffi- 
ciently small (for details, see Fig. 2). A sketch as how to 
find the conduction eigenchannels in this case has been 
given in Sec. IV. E. Below we complete that discussion by 
presenting a mathematical justification of the method. 

We proceed in three steps: (i) since a position of Sl 
and Sr planes in Eq. ([9]) does not meet an asymptotic 
limit criterion, we expand the Green's function in terms 
of scattering Bloch states in contrary to expansion (fTl| 
involving unperturbed states; (ii) we use the Kubo for- 
mula ([9]) and express conductance as a trace of the appro- 
priately defined current operator in fc-space; the obtained 
operator is identified with transmission matrix t't; (iii) 
finally, we build up the equivalent site-angular momen- 
tum representation of the transmission matrix, that leads 
us to the required formulae [Eqs. ([29"jl. ([30]) of Sec. IV.E]. 



1. Expansion of the Green's function. 

To begin the proof, we note that the asymptotic repre- 
sentation ([T4l) of the Green's function can be rewritten as 



G+(r, r',£:) | Z ^ TOO (Al) 



kk' 



where an expansion is performed over perturbed Bloch 
waves $ fc (r,£:) = $ k (r,E)/y/Vk~ and $ k (r,E) = 
^ K (r , E) / y/vk which were introduced in Sec. III. We 
remind, that function 3>fc(r, E) is the solution of the 
Lippmann-Schwinger equation ^ associated with the 
in-coming unperturbed Bloch wave ^(r, E) / ^/vk in the 
left lead propagating towards atomic constriction (vj- oc 
v k > 0), while the second function, <I>fc(r, E), is the so- 
lution of the equation associated with the in-coming un- 
perturbed Bloch wave (r, E)j y/V-k in the right lead 
{v- k oc v z _ k < 0). 

Consider now a general case, when z = zl and z' = zr 
points related to planes Sl and Sr (see Fig. 2) are not 
taken infinitely far from atomic constriction. We note, 
that for the conductance evaluation one needs only the 
back-scattering term 5G + (r, r', E) in Eq. |T]), which is 
the solution of the Schrodinger equation without a source 
term. Therefore, we can expand SG + (r, r', E) over the 
eigenf unctions of the whole system corresponding to en- 
ergy E. These are the propagating perturbed Bloch 
states, 



and the evanescent states, 

X a (r,E)=J2r a , nL (E)ct>L(r--R n ,E), (A3) 

ii L 

where, in turn, both type of functions are expanded over 
atomic orbitals, and Ck,nL and T a ^„L are expansion co- 
efficients. 

In particular, matrix Q K ^ n L is related to the matrix 
C Ki nL corresponding to the unperturbed Bloch waves (for 
details, see 

r _-| -i nn' 

e k ,nL(E) = [{1 - Go(E) At(E)} \ LL C k , n , L >{E) 

n' L' 

where the Ai-matrix describes the whole scattering re- 
gion (a vacuum barrier plus an atomic constriction) 
and Go is the structural Green's function of the three- 
dimensional periodic bulk crystal. 

The evanescent states can be, for example, the surface 
states which are perturbed by an atomic constriction, as 
well as the states localized around impurities. 

Both equations (|A2|) and (|A3[) can be joined to one 
matrix equation: 



(J,E), 



where (...) stands for a column- vector, while [. . . ] de- 
notes a matrix. We also expand the conjugated states as 



( *(r,£)\ _ 


' 6(25) " 


\X(v,E) ) 





*(r,£) 
X*(r,E) 



)- 


' 6(2?) " 




T(E)_ 



(r,E). 



We expand the back-scattering term 6G + (r, r', E) over 
eigenfunctions ^(r, E) and X a (r,E) with energy E: 

5G + (r,r',E)\ 
= E £<^(r,£) GT L \{E) E) 

neS L ,n'eS R LL' 

' " F 00 (E) F 01 (E) 
F W (E) F n (E) 

J2*k(r,E)F™(E)* k ,(r',E) 

kk' 

J2*k(r,E)F$(E)X f3 (r',E) 

fc/3 

J2x*(r,E)F™,(E)* k ,(r',E) 

ah' 

J2x*(r,E)F$ ) (E)X p (T,E) 



*(t,E) 
X*(r,E) 



V{r',E) 
X(r',E) 



= SG^\r,r',E)+dG^(r,r',E). 



(A4) 



$ k (r,E) = Y, Q k,nL(E) 0i(r-R„,S), 



(A2) 



Only the first term 
6G^(r,r',E) = £ * k (r,E) F™,(E) * k ,(r>, E) (A5) 



kk' 
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exists in the asymptotic case, while other three terms, 
denoted as 6G\ (r, r , E), contain contributions of the 
evanescent states which decay within the leads. 

The unknown matrix F — {i™} in Eq. (|A4[) satisfies 
the equation (below we will skip for simplicity the obvious 
energy dependence): 



F oo 


F oi 




e 


pio 


F n 


X 




r 



Glr, 



where symbol T denotes transpose, and matrix Glr is 
introduced which contains selected matrix elements of 
the structural Green's function G™ l L , with n € Sl and 
n'eS R . 

To find matrices F lJ , we introduce the block-matrix 
[fit At] which is a solution of the following equation 
(here symbol f denotes Hcrmitian conjugate): 



est 


CAt 




1 


T0t 


rAt 




1 



(A6) 



This means, that matrix \B' A*] is the inverse or (in 
general case) the pseudoinverse^ matrix to the matrix 
with 6- and T-blocks^. 

The equation similar to Eq. (|A6[) exists for other ma- 
trices: 



B* 
A !i 



qT t t 



B*e T b* r T 
A*e T A*f T 



1 o 
o 1 



(A7) 

where symbol * denotes complex conjugate of the matrix 
elements only. 

With the use of Eqs. (|A6|) . (IA7p we find solution for 
F = {FV}: 



' F 0Q 


poi ' 




' B*G L R0 t 


B* G LR At " 


plQ 


F u 




_ A*G LR 23+ 


A*G LR At 



Then expansion coefficients of the SG^ term read as 

= E T,^n L (E)[GME)]ZK'L,A E )- 

n^S-^.n 1 G5r LL' 



With help of matrix r L p (E), 



2m 



(A8) 



(with positive velocities under square-roots) we write 
down Eq. (|A5|) as 



6G { ^(r,r',E)\ 
= -2m^ h (r,E) [r^(E)] kh , $ k >(r',E). 



(A9) 



Comparing obtained equation with one for the asymp- 
totic limit [Eq. ([Aip ] and taking into account that in 
fact r£^ (E) is independent on the positions of the L and 
R planes, we obtain: 



[^(^)] fefe '=[-- 1 (^)] fefe ' 



Evaluation of conductance 



(A10) 



When conductance is evaluated, only 5G+(r,r',E) 
term in expansion (j A4|) contributes. Other term, 



SG + (r,r',E), contains the evanescent states which de- 
cay within the leads and, therefore, have zero velocities 
along the current flow. When Eq. (|A9[) is inserted into 
the Kubo formula J9l) we obtain: 



9 = 9o E E 

where all matrices are evaluated at the Fermi energy 
(Ep = E k = E k i = . . .). Here we have introduced the 
current operators related to the left (L) lead, 

= 2tt J dS 
and to the right (R) one 

'kk 



1 



tlr 



kk' 



^k'ki 



1 

T 



(All) 



$ k *(r,E)id z **{t,E) 



Jkk'(E) 



= 2ir J dS 

Si, 



$> k (r,E)id z 9fr{T,E) 



(A12) 
(A13) 



where basis v of eigenchannels (see Sec. Ill) was used to 
evaluate the matrix elements. With help of (|A12[) and 
(|A13[) . equation (|A11[) for conductance g takes a form: 

9 = g Tr fc (rtrO LR ) (AH) 

- SoTr fc (tV) = g Tr k (j L ) , 

8 kk > is the unitary operator in the fc-space: 
-l 



where = 



O 



LR 



(A15) 



1. 



and Eqs. (|A10j) and (|A13|) were used. 

When the perturbed Bloch waves are expanded over 
atomic orbitals, 



*h(r,£?) 



1 



2vr 
1 



v k 
1 



kk' 



*fe(r,£) 

J2 ^2x k!nL (E)<Pl(r,E), 

n£S R L 

^ ^X k , nL (E)r L (r,E), 



2lT u&Sl L 
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current operators J 1 " and J7 R take a form (energy de- 
pendence is skipped): 



= e e*um^:w> ( Aie ) 

uu'gSl LL' 

^kk' = E E Xfe -" L ^lI' K'L'k" 



nn'^Sft. LL' 

where operators X and X are defined as 



2tt 



1/2 



iMi(^) = ( — ) e k , nL (E), (A17) 



1/2 



Sfe,ni(-2')) 



and operators Z?l and Dr were introduced in Eq. (|2T 



3. Decomposition of current operator 



According to Eq. (|A14|) conductance is determined by 
operator J 1 " = t^t. The formal representation (|A16[) 
of the current operator J 1 " in fc-space is not suitable 
for practical implementation. It can be used, however, 
to find an equivalent representation in (n, L)-space. Us- 
ing (|A16|) together with representation of operator 
[Eqs. (23), (24)] we can decompose J7 L in two terms: 



J 



(A18) 



where 



= +X* U L 



A+ 




U[[X* 



(A19) 



has only positive (non- negative) eigenvalues, while 




ul[x*] 



(A20) 



has only negative (or zero-) eigenvalues, and are diag- 
onal matrices with nonnegative elements [Eqs. (23), (24)]. 

The advantage of operators over operator Jx, is that 
matrices are positive definite, while anti-symmetric 
matrix is not. As we show below in Appendix B, that 
makes possible to find for operators an equivalent, 

site-angular momentum (n, L) representation, while it is 
not the case for full operator We mention also, that 
in the asymptotic limit the contribution due to ' van- 
ishes. However, in general case, spectra of both operators 
should be found. 

To clarify these ideas, let us elucidate a physical mean- 
ing of Eq. (|A18[) . Consider new basis functions in (n, L)- 
space which are constructed with help of unitary ma- 
trix J7l, 



Xa (r,E)= E^(r^) [UL(E)] n 



L,a> 



so that the velocity operator Z?l = UlD^U^ in {n,L)- 
space introduced in Eq. (21) is diagonal in the new basis: 



[DUE)} 



a/3 



(A21) 



dS 



X* a {T,E)id zX ^,E) 



Let us examine further the wave function < &fc(r, E) which 
defines the current operator [Eq. (|A12[) ]. This func- 
tion is the solution of the Lippmann-Schwinger equation 
associated with the in-coming unperturbed Bloch wave 
Q^*(r, E) (initial channel) propagating from the right 
(R) lead towards the nanocontact and being scattered 
on it. Within the left (L) lead (r e Sl) we have: 



$ fc (r,£) = ^ k + \v,E) + $ { k \r,E), 



where 



^2xa(r,E)(X*U h )l k , (A22) 



here the (a ± )-sums run over the half of the indices of 
the basis functions Xa( r iE), corresponding either to the 
"positive" or to the "negative" window of the spectrum. 
Thus, the wave function < &fc(r, E) in the L lead is a linear 

combination of two functions. The first one, <I> fc + (r, E) 
carries the flux in the initial direction of the in-coming 
wave with momentum — fc, from the L electrode to the R 

one. The second function, & k \r,E), carries the flux in 
the direction being opposite to the propagation direction 
of the in-coming wave. An example illustrating this idea 
for free electrons is shown in Fig. 3. 

One can check, that operator jS is defined on func- 



tions $> k ' (r, E) only, while ' is defined on & K ' (r, E): 



[^t\ E )] kk '=^ / dS 



$ k ±) ~(r,E)id z W« J (r,E) 



because, according to Eqs. (|A21[) and (|A22[) . the 



cross-terms involving both functions, ^ k + \v, E) and 

\r, E), vanish. 
Obviously, that if the 5l plane is chosen way behind 
the scattering region, the & k \r,E) turns to zero and 



the contribution due to 
both operators, 



(-) 



vanishes. In general case, 



and jS , should be considered. 



n£S L 



4. Spectrum of current operators. 

In this section we show that transmission probability of 
the conduction eigenchannel v is given by T v = r+ — r~ , 
where r+ and (— t~) are positive and negative eigenval- 
ues of operators and , respectively. 
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Assume, that we know matrix 9 (introduced in Sec. Ill) 
which is the unitary transform to the basis of eigenchan- 
nels. According to Ref. [4?] we have: 

¥„(r,25) = J2^(r,E)9UE) 



= $ 



k 

(+) 



(r,E) + $i \r,E). 



Thus, we have decomposed the wave function of each 

eigenchannel in two terms: (fr^ = J2k @kv 

Consider further the transmission matrix T in the basis 
of eigenchannels: 



2ir / dS 

t(+) + t(") 



with 



T<±) = 2tt / dS 



>k) *(r,E)id z <S>f\r,E) 



Operators are related to current operators via 
the unitary transform 9: = 9j£ ±) 9'<. 

The unitary transform 9 diagonalizes the full operator 
J L = J"l +) + Jl ^- Using the representation (|AL6[1 of 
operators we can check, that jj^' and do 

not commutate. Thus, the T^' operators do not have a 
diagonal form in the basis of eigenchannels. Therefore, 
we represent each of these operators as sum of diagonal 
and off-diagonal terms: 



T (-) 



T 



.(-) 
(+) 



J T (-) 

I'll ! 



(A23) 



here 8T„„ ' = 0, and off-diagonal contributions are of 
different signs because the sum of two matrices is the 
diagonal matrix T U S^. The diagonal terms are 



»(+) 



T+S 



□ (-) 
T 



(A24) 



where 



= 2vr / dS* 



and transmission probability of the eigenchannel is 
given by T v — r+ — t~ . Here the first term arises due 
to all multiple scattering contributions in the direction 
of the current, while the second term is due to scattering 
contributions in the opposite direction. 

In spite of the fact that operators T^ + ) and T^~> can 
not be diagonalized within the unique unitary transform, 
one can show that off-diagonal terms ST^ ' ^ are small. 



They are determined by functions (r, E) with nega- 
tive velocities with respect to the initial velocity of the 
in-coming channel v incident from the right (R) lead. The 



function <f> 



(-) 



(r, E) in the left (L) lead is small as com- 

<+)/ 



pared with & v (r, E). When the Sl plane is placed far 

enough from the surface, the "* (r, E) collects all mul- 
tiple scattering events in the opposite direction to the 
current flow. However, such scattering processes are pos- 
sible only due to small inhomogeneities of the potential 
which is not exactly the bulk one around zl (see Fig. 3). 
Thus, we can write down: 



\§l\r,E)\ 



p\^l + \r,E)\ 



where parameter p < 1 has a meaning of reflection am- 
plitude (see caption of Fig. 3). The p 2 is of the order 
of reflection probability, which (in realistic calculations) 

is p 2 ~ t~ 7t+ ~ 0.1. Thus, off-diagonal terms, ST^\ 
contain small parameter A ~ p 2 ~ 0.1. Since 5Tj^/ = 0, 
according to the perturbation theory^ the difference be- 

m . o(±) 

tween eigenvalues of operators Tp, v ' [Eq. (|A23| )] and T vfJL 

[Ea. (|A24p ] appears only in the second order, which is 
0(A 2 ) ~ 0.01. Therefore, the operators provide 

a feasible way to find Tu with good enough precision. 

To find the spectrum of positive definite operators 
c7l we represent them in the form £j_E±, where £± 
maps the site-orbital space on the fc-space. We can prove 
(see Appendix B) that spectrum of £lj_£± is the same as 

for operator E±£!|_ which is defined in (n, L)-space and 
reads as 



T(±) 



1/2 



[X*]*O hR X* (±D 



1/2 



Here we have used the unit operator in fc-space C LR = 

t lr ( t rl)_ = ^fefe' introduced in Eq. (|A15|) . By using 
equations (|A8|) . (|A15|) . (|A16|) . (|A17|) . we obtain expres- 
sion for operator T^': 

2"(±) 



± ±D 



±\l/2 



[X* 



± l±fr 



.1/2 



BG 



± ±D 



t 

RL 

t\V2 



e 

[ir] t ) e* (±d 



/ t 



X* (±D 



1/2 



LR 
±\l/2 



Glr^rG^ (±D 



1/2 



where orthogonality relation^ B^C = [B*]^C* — 1 was 
used. Thus, we come to Eq. (30) of Sec. VI. 



APPENDIX B 

Let us consider operator A fcfc / (E) acting in fc-space of 
a following kind: 

A = XAXK 
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where A is a positive definite hermitian operator acting 
in the site-orbital (n, L)-space. Let Q v k be a unitary 
transform which diagonalizes operator A, 

6A<d f = diag{A„} = A . 

Our goal is to find a spectrum of A with help of (n, L)- 
space. Following Sec. Ill, we represent A in the form: 

A = S^E = XAX\ 

with S = A 1 / 2 ATt and £t = ATA 1 / 2 . A square-root from 
the positive definite hermitian operator is defined as 

A l/2 = (A l/2 )t = UA l/2 L/t; 

where the unitary matrix U transforms A to the diagonal 
form A = W AU. 

Let us show that a spectrum of operator A = S^S act- 
ing in the fc-space, is the same as a spectrum of operator 
SS 1 * acting in the (n, L)-space. Let be such matrix 
that 

= fiSet = diag{V^} = V 
Solution for £1 is £l v , nL = J2k diag{l/VA7}Oi,fc£i. nL . 



Matrix £1 diagonalizes EEL 

= v 1 (es t se t )(es t se t )v 1 

= A^^oAoV 1 = diag{A„}. 

To finish the proof, one has to check that £7 is indeed 
the unitary transform, i.e. the following properties should 
hold: (i) fifit = and (ii) fl^fl = & nn '&LL'- First, we 
check property (i) : 

= A^esW)^ 1 = Ao'AoAq^I. 
Next, we check property (ii): 

= S(e t Ao 1 6)S t = £[EtE] _1 £t = 1, 

so the proof is complete. 
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and E„ e s L Si r a ,„LAf i (3 = 5 a p. In order the n-sum 
be converged, matrix elements of C, B\ F and A" 1 " should 
carry coefficients ~ \f\j N s (where N s = N x N y is num- 
ber of atoms per cross-section of the Born-von Karman 
supercell). On the other hand, since matrix A 1 "] is the 
inverse to the one with C and Y blocks, the following rela- 
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The second term in this equation is zero, because the num- 
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~ N x N y = N s , so that J2 a Ai iiQ r ain , L , ~ N a /N s -> 0. 
Finally, we obtain: ^ fe B\ L k &k,nL = 5 nn >S LL >. 



